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We analyze the propagation of elastic waves in soft materials subjected to finite deformations. We derive 
explicit dispersion relations, and apply these results to study elastic wave propagation in (i) nearly incompress¬ 
ible materials such as biological tissues and polymers, and (ii) negative Poisson’s ratio or auxetic materials. 
We find that for nearly incompressible materials transverse wave velocities exhibit strong dependence on 
direction of propagation and initial strain state, whereas the longitudinal component is not affected signif¬ 
icantly until extreme levels of deformations are attained. For highly compressible materials, we show that 
both pressure and shear wave velocities depend strongly on initial deformation and direction of propagation. 
When compression is applied, longitudinal wave velocity decreases in positive bulk modulus materials, and 
increases for negative bulk modulus materials; this is regardless the direction of wave prorogation. We demon¬ 
strate that finite deformations influence elastic wave propagation through combinations of induced effective 
compressibility and stiffness. 
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The propagation of elastic waves has been investigated 
intensively^“^^ because the understanding of the phe¬ 
nomenon is vital for a large variety of applications from 
non-invasive materials testing and medical imaging for 
health care to petroleum exploration. Recently, the field 
of acoustic or phononic metamaterials has attracted a 
considerable attention. The peculiarity of these meta¬ 
materials originates in their microstructure, which can 
be tailored to give rise to various effects such as local 
resonances^^, band-gaps^^ and cloaking^^. Furthermore, 
soft metamaterials, due to their capability to sustain 
large deformations, open promising opportunities of ma¬ 
nipulating acoustic characteristics via deformation^^“^^. 

In this work, we derive explicit dispersion relations 
for finitely deformed materials by means of the su¬ 
perimposed on large deformations Bloch wave anal¬ 
ysis frequently used for investigation of elastic wave 
propagation^’^^’^^’^^’^^. The availability of the explicit 
dispersion relations is important for designing meehan- 
otunable acoustic metamaterials. Moreover, the infor¬ 
mation may greatly benefit non-invasive medical diag¬ 
nostic techniques by providing the important informa¬ 
tion on dependence of elastic wave propagation on pre- 
stress/pre-strain conditions, which is a common state 
of biological tissues. By application of the derived ex¬ 
plicit expressions, we show the role of deformation in 
wave propagation in soft media undergoing finite defor¬ 
mations. Moreover, we extent the analysis to a class of 
exotic metamaterials characterize by negative Poisson’s 
ratio (NPR) behaviour. Examples of the NPR materi¬ 
als, also known as auxeties, include living bone tissue^^. 
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skin^^, blood vessels^^, certain rocks and minerals^^, and 
artificial materials^^. As we shall show, elastic wave 
propagation in these materials is significantly affected by 
deformation. 

To analyse the finitely deformed state, we introduce 
the deformation gradient F(X, t) = Vx (8) x(X, t), where 
X and X are position vectors in the reference and current 
configurations, respectively. To take into account finite 
deformation non-linear effects as well as material non¬ 
linearity, we analyse the wave propagation in terms of 
infinitesimal plane waves superimposed on a finitely de¬ 
formed state^’^. To account for the stiffening effects (due 
to, for example, finite extensibility of polymer chains, 
or due to collective straightening of collagen fibers in 
biological tissues) in finitely deformed media, we make 
use of the strain-energy density function corresponding 
to an approximation of Arruda-Boyce model^^, namely 
the Gent model^^ 
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where /i is the initial shear modulus, K is the bulk mod¬ 
ulus, /i = tr B is the first invariant of the left Cauchy- 
Green tensor B = FF^ and J = det F. The model neatly 
covers the stiffening of the material with the deforma¬ 
tion; as the deformation attends the level of /i = 3 -h 
the energy function becomes unbounded and a dramatic 
increase in stress occurs. Consequently, Jm is locking 
parameter. Clearly, when ^ oc, the strain-energy 
function (1) reduces to the neo-Hookean one, namely 

^(F) = |(/i - 3) - Min J + {KI2 - m/3) (J - l)^. (2) 
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Acoustic tensor for finitely deformed neo-Hookean ma¬ 
terial ( 2 ) takes the form 

Q(n) = ain (g) n -h a 2 (I — n (g) n), (3) 

where the unit vector n defines the direction of propa¬ 
gation of the wave; n (g) n is the projection on the di¬ 
rection n; (I — n (g) n) is the projection on the plane 
normal to n; ai = (iC — 2/i/3)J + ^(1 + n-B-n) and 

a 2 = (n-B-n). Consequently, there always exist one 

longitudinal (pressure) and two transverse (shear) waves 
for any direction of propagation n. Phase velocities of 
these waves can be calculated as 

Cpw = 'Jai/p and = sfa^P-, (4) 

where p is the material density in the deformed state. In 
the limit of small deformation s the results re duce to the 
well known relations Cp^ = and Csw = 

\fpTp- 

We present the results in the generic form where possi¬ 
ble, however, the illustrative examples are given for uni¬ 
axial tension 


F = Aei (g) ei -h A(e 2 (g) 62 + 63 (g) es), (5) 


where A is the applied stretch ratio and A = 
is defined through A and compressibility of the material. 
Note that the compressibility of the material is defined 
by ratio Kjp. In the limit of linear elasticity the elastic 
moduli are related through 


K _ 

TT “ 3(1 -2z/)’ 


(6) 


where v is Poisson’s ratio. Thus, —1/3 < K/p < 00 with 
/i > 0. Note that for —1/3 < K/p < 0 the material is 
stable only if constrained^^“^^ . Besides, it follows from 
(6) that matter exhibits auxetic behaviour when —1/3 < 
K/p < 2/3. 

For nearly incompressible materials {K/p ^ 1 and A 
A“^/^) the dependence of the pressure wave velocity on 
the direction of propagation and initial stress state is 
relatively weak. However, the velocities of shear waves 
depend strongly on the direction of propagation and the 
deformation state (FIG. 1). In particular, we observe 
prominent minima and maxima of the velocities when 
waves propagate along the principal directions. Here and 
thereafter the velocities are normalized by the value in 
the undeformed state. 

In contrast to nearly incompressible materials, for 
highly compressible matter, the velocity of pressure wave 
depends strongly on the direction of propagation and 
initial stress state (see FIG. 2 (a)). More specifically 
the velocity of pressure wave increases when material is 
stretched and decreases when it is compressed. Remark¬ 
ably, if K/p varies in the vicinity of 2/3 which corre¬ 
sponds to Poisson’s ratio u in the vicinity of 0, then the 
velocities of pressure and shear waves are not influenced 



FIG. 1. Polar diagram of the phase velocities for the shear waves 
{K/fi » 1 ). 


by the initial stress-strain state when wave vector n lies 
in the plane of orthotropy orthogonal to the axis of ten¬ 
sion (ei), see FIG. 2. Propagation of shear waves in 
highly compressible media differs significantly from the 
one in nearly incompressible materials. In particular, un¬ 
der compression the velocities of shear waves decrease in 
any direction if K/p < 1 (FIG. 2 (b)), while it can ei¬ 
ther decrease or increase depending on the propagation 
direction in the nearly incompressible case (see FIG. 1). 



FIG. 2. Polar diagrams of the phase velocities for the pressure (a) 
and shear (b) waves, K/fi = 0.7. 


Next, we consider the auxetic materials with the nega¬ 
tive bulk modulus. Yet can such materials exist? Would 
they be stable? How to construct them? These ques¬ 
tions have recently arisen in many papers^^’^^’^^“^^ and 
still the topic is open for discussion. Wang and Lakes 
in their article^^ report that bulk modulus can be varied 
within the range — ^ < iF < oc. Based on this and our 
previous estimations, we examine the material behaviour 
when —1/3 < K/p < 0. An example for an auxetic ma¬ 
terial with K/p = —0.3 is shown in FIG. 3 for pressure 
(a) and shear phase velocities (b). We observe that the 
velocity of the pressure wave increases in any direction 
of propagation n when the material undergoes compres¬ 
sion (FIG. 3 (a)) in contrast to the case of positive bulk 
modulus (see FIG. 2 (a)). Moreover, the velocity of the 
pressure wave increases significantly and has maximum 
when wave vector n lies in the plane orthogonal to the 
axis of compression. Secondly, velocities of shear waves 
(FIG. 3 (b)) decrease considerably for all directions of 
propagations n in contrast to the effect of deformation 
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observed in the nearly incompressible materials (FIG. 1). 




FIG. 4. Polar diagram of energy curves for the shear waves 
(K/m » 1). 


FIG. 3. Polar diagrams of the phase velocities for the pressure (a) 

and shear (b) waves, Kjii = —0.3. form 


The explicit expressions for corresponding angular fre¬ 
quencies are 


U). 
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( 7 ) 

(8) 


where k is wave vector, k = |k| is wave number and po 
is the density of the undeformed material. Note that 
the wave vector can be also mapped into the undeformed 
configuration as k^ = J“^F^k. 

The transmission velocity of a wave packet or the group 
velocity can be found by the use of the formula Vg = 
VkCo’ O 6uj = dli ' Vg. Thus, expressions (7) and (8) 
yield 

Mk + B.k) + (X-f)J^k 


and 


To illustrate expressions (9) and (10) we plot the en¬ 
ergy curves in FIG. 4. Usually for homogeneous isotropic 
materials the velocity and energy curves coincide^^’^^. 
However, in case of the nearly incompressible pre-stressed 
homogeneous isotropic material we find that velocities of 
energy and waves propagation coincide only for principal 
directions and when wave vector n lies in the plane or¬ 
thogonal to the axis of tension. For other directions of 
propagation n they are significantly different (compare 
FIG. 4 and FIG. 1). 

Next we examine the influence of deformation on wave 
propagation in materials characterized by strong stiffen¬ 
ing effects. In particular, we make use of the Gent ma¬ 
terial model (1); the corresponding acoustic tensor takes 


Q(n) 




+ (^(n • B • n)I + 2n 0 n : B 0 , 


( 11 ) 


where ^ = 3 + — U and B (g) is isomer^^ of 

tensor of fourth order, i.e. di (g) d 2 (g) ds (g) = 

di (g) dj (g) d/c (g) ds, where (ijks) is a permutation of the set 
(1234). One can see from (11) that in general case the 
waves are neither pure longitudinal nor pure transverse in 
contrast to the waves prorogating in neo-Hookean mate¬ 
rials. However, if wave vector n lies in one of the planes 
of orthotropy then we always have one pure transverse 
(shear) wave with the velocity 


^sw 


\J ^2Jml (pO • 


( 12 ) 


Note that expression (12) is general and valid for any 
finite deformation F and compressibility Kjp. 

Similarly to the neo-Hookean results, for waves prop¬ 
agating along principal directions in Gent material, we 
have only pure modes. The corresponding expressions 
of velocities can be derived directly from (11), however 
they are rather tedious. For waves propagating along 
the direction of loading (n = ei) in nearly incompress¬ 
ible matter the velocities of longitudinal (pressure) and 
transverse (shear) waves reduce to 




K 


\3 Jrr, 


+ 


JmA(A3 + (3+J„)A-2)^ 


(13) 


and Cg^ju — Jmk'/ 

It is easy to see from (13) that for nearly incompress¬ 
ible materials the dependence of the longitudinal wave 
velocity is relatively weak unless extreme levels of defor¬ 
mations are attained. At these extreme deformations, 
the stiffening effect manifests in a significant increase of 
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the effective shear modulus, which becomes comparable 
with the bulk modulus. We illustrate this dependence of 
the longitudinal wave velocity on the direction of propa¬ 
gation in FIG. 5. We observe a dramatic change in the 
velocity profile of the longitudinal wave in extremely de¬ 
formed materials. Moreover, the velocity of longitudinal 
wave increases in both compressed and stretched materi¬ 
als. 



FIG. 5. Polar diagram of the phase velocity for the quasi-pressure 
waves in Gent material (Jm = 3, Kj\i = 300). 

For the transverse waves, we find that the velocities 
depend strongly on the deformation and the stiffening 
parameter (14). Examples of the dependence of the shear 
waves propagation on deformation for Gent material with 
= 3 and Kj are shown on FIG. 6. The results 
clearly indicate the significant effect of the deformation 
induced stiffening (compare FIG. 1 and FIG. 6). We 
observe that the dependence of the shear wave velocity 
on deformation and prorogation direction increases when 
the locking parameter decreases (which corresponds 
to earlier stiffening of the material with deformation). 
Gomparing the quasi shear wave velocity profiles with 
the pure shear ones, we observe that the velocity of quasi¬ 
shear wave has maxima for the non-principal directions 
(see FIG. 6 (a)). 



FIG. 6. Polar diagrams of the phase velocities for the quasi-shear 
(a) and pure shear (b) waves in Gent material (Jm = 3, Kj\i ^ 1). 

We derived the explicit dispersion relations for finitely 
deformed materials, and demonstrated the significant 
role of the deformation on the elastic wave prorogation on 
examples of both nearly incompressible and extreme aux- 
etic materials. These findings may guide further design 


of mechanotunable acoustic metamaterials and phononic 
crystals with a large range of constituent properties. The 
local distribution of strain field in these engineered ma¬ 
terials can be used to induce regions with extreme varia¬ 
tion of the local phononic properties to give rise to var¬ 
ious acoustic effects. This opens a very rich and broad 
research avenue for designing tunable acoustic/phononic 
metamaterials. 
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